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Abstract. Motivated by recent scanning tunneling and photocmission spectroscopy 
measurements on self-organized gold chains on a germanium surface we reinvestigate 
the local single-particle spectral properties of Luttinger liquids. In the first part we 
use the bosonization approach to exactly compute the local spectral function of a 
simplified field theoretical low-energy model and take a closer look at scaling properties 
as a function of the ratio of energy and temperature. Translational invariant Luttinger 
liquids as well as those with an open boundary (cut chain geometry) are considered. We 
explicitly show that the scaling functions of both setups have the same analytic form. 
The scaling behavior suggests a variety of consistency checks which can be performed 
on measured data to experimentally verify Luttinger liquid behavior. In a second 
part we approximately compute the local spectral function of a microscopic lattice 
model — the extended Hubbard model — close to an open boundary using the functional 
renormalization group. We show that as a function of energy and temperature it 
follows the field theoretical prediction in the low-energy regime and point out the 
importance of nonuniversal energy scales inherent to any microscopic model. The 
spatial dependence of this spectral function is characterized by oscillatory behavior 
and an envelope function which follows a power law both in accordance with the field 
theoretical continuum model. Interestingly, for the lattice model we find a phase 
shift which is proportional to the two-particle interaction and not accounted for in 
the standard bosonization approach to Luttinger liquids with an open boundary. We 
briefly comment on the effects of several one-dimensional branches cutting the Fermi 
energy and Rashba spin-orbit interaction. 
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1. Introduction 

Over decades theoretical studies of the single-particle spectral properties of metallic 
one-dimensional (Id) correlated electron systems — so-called Luttinger liquids (LLs)— 
were ahead of the experimental attempts to find or synthesize appropriate quasi Id 
materials and perform spectroscopy on them. In fact, while at the begining of the 
1990's a clear picture of the basic spectroscopic properties of translational invariant LLs 
was established (for reviews see e.g. Refs. PQ, [2] and [3]) this period witnessed the first 
serious attempts to experimentally verify the specific spectroscopic signatures of LLs 
[3]. These are the (i) low-energy power-law suppression of the local spectral function 
p{uj) ~ \uj\ a for energies u close to the chemical potential El [7] with a depending 
on the two-particle interaction as well as (ii) the appearance of two dispersing features 
in the momentum resolved spectral function p(k, to) (spin-charge separation) [HI El \10\ 
instead of a single quasi-particle peak of a Fermi liquid. For finite temperatures T the 
suppression of the spectral weight as a function of u is cut off by T and one finds the 
scaling behavior p ~ T a S a (u/T) with a a-dependent scaling function S a in which the 
two energy scales u and T only enter via their ratio [TT] . These results were exclusively 
obtained using bosonization within the Tomonaga-Luttinger (TL) model [U El E] • 

Using the modern language of renormalization group (RG) methods the 
(translational invariant) TL model is the exactly solvable effective low-energy fixed 
point model for a large class of metallic Id correlated electron systems — the LLs [12j . 
It thus plays the same role as the free Fermi gas in Fermi liquid theory. The model 
has two strictly linear branches of right- and left-moving fermions and two-particle 
scattering is restricted to processes with small momentum transfer \q\ <C kp, with the 
Fermi momentum kp. These processes as well as the kinetic energy can be written as 
quadratic forms of the densities of right- and left-moving fermions which obey bosonic 
commutation relations. In most calculations in addition the momentum dependence of 
the low-momentum scattering processes #2 and g± are neglected and momentum integrals 
are regularized in the ultraviolet introducing a cutoff 'by hand' (for an exception see 
Ref. [9]). One can extend the resulting scale- free, field theoretical model by allowing 
for additional two-particle scattering processes. These turn out to be RG irrelevant in 
a wide parameter regime [13]. The most important of these processes is the so-called 
gi j_-process (in the g-ology classification [13]) with momentum transfer 2fcp between 
two scattering fermions of opposite spin. 

In some of the early experiments on Id chains these were obviously interrupted 
by local impurities [4]. A simple model of an inhomogeneous LL is the open boundary 
analog of the TL model. Interestingly, a LL is very susceptible towards single-particle 
perturbations with momentum transfer 2kp [7J and on asymptotically low energy scales 
even a single weak impurity has the same effects on the spectral properties as an open 
boundary [T3]. Triggered by this theoretical insight and the early experiments, the 
spectral properties of the open boundary analog of the TL model were studied [TBI [16] . 
The local spectral function close to the boundary shows power-law behavior as a function 
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of u but with an exponent «b different from the bulk one a. As in the translational 
invariant case in this model only those low-energy scattering terms are kept which can 
be written as quadratic forms in bosonic densities. Only recently it was shown that 
a large class of further two-particle processes appearing in a Id system with an open 
boundary are indeed RG irrelevant [T7] . 

The latest scanning tunneling spectroscopy (STS) and photoemission spectroscopy 
(PES) measurements on different classes of Id metallic systems [T8l [T9| [20| [2T| [22] 
impressively demonstrated that the experiments caught up and more refined questions 
must now be answered by theory. Important ones are: How do effects which are not 
captured by the low-energy fixed point model, such as the momentum dependence of the 
two-particle interaction and the nonlinearity of the single-particle dispersion influence 
the spectral functions? What is the energy scale of a given microscopic model on which 
the low-energy LL physics sets in? How do scaling functions for lattice models look 
like in detail? Here we shed some light on the last two questions and briefly comment 
on the first one. It is widely believed that neglecting the momentum dependence of 
the interaction and regularizing momentum integrals in the ultraviolet 'by hand' has 
no effect on the low-energy physics of LL's. This is indeed correct if all energy scales 
are sent to zero, that is for p(u) and p(±k-p,uj): at small u the spectral properties are 
unaffected by the details of the momentum dependence of the g's. However, if p(k, u) as 
a function of u is studied at fixed k 7^ ±kp, as it is usually done in angular resolved PES, 
details of the momentum dependence of the interaction matter. This was investigated in 
Ref . [23] . An overview on the effects of the nonlinearity of the single-particle dispersion 
can be found in the very recent review Ref. [23]. 

This paper is organized as follows. In Sect. [2] we compute the local spectral 
function of the translationally invariant and the open boundary continuum TL model 
using bosonization. We show that both display scaling in u/T and that the scaling 
functions have the same analytic form. We next compute the spectral function of 
the extended Hubbard model on the lattice close to an open boundary as a function 
of energy, temperature and position in Sect. [3] For this an approximate method is 
used which is based on the functional RG approach [25]. It is devised for weak to 
intermediate two-particle interactions. In particular, we concentrate on inhomogeneous 
LLs as the boundary exponent «b characterizing the spectral function close to an open 
boundary is linear in the two-particle interaction while the bulk exponent is quadratic 
(see below). Varying the microscopic parameters of the extended Hubbard model we 
can tune the strength of the different scattering processes and thus study the crossover 
between nonuniversal behavior and the low-energy LL physics. We perform a scaling 
analysis of the spectral function as a function of u and T and show that the spectral 
weight close to the boundary follows the bosonization prediction within the universal 
low-energy regime. The position dependence of the spectral function is characterized 
by oscillatory behavior and a power-law envelope function in accordance with the result 
for the TL. Interestingly, we additionally find a phase shift which is proportional to the 
two-particle interaction and not accounted for in the standard bosonization procedure. 
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We summarize our results in Sect. H]and briefly comment on how spin-orbit interaction 
and several bands crossing the Fermi surface — both being potentially important effects 
in recent experiments — influence the single-particle spectral functions. 



2. Scaling functions of the Tomonaga-Luttinger model 

In this section we derive closed analytic expressions for the single-particle Green function 
and the related local spectral function of the TL model with and without an open 
boundary at finite T. We then closely inspect the scaling form of the spectral functions. 

In field theoretical notation (see e.g. Ref. [26]) the Hamiltonian density of the TL 
model in spin-charge separated form reads 



16tt 



-^(d x ^) 2 + K v (d x Q v ) 2 



(1) 



with the canonical Bose fields and their dual fields Q u . Within the TL model the 
charge and spin velocities f CjS as well as the LL parameters K CjS are free parameters. 
If the model is used to describe the low-energy physics of an underlying microscopic 
model they become functions of the corresponding model parameters and the band 
filling [H El E] • For spin- rotational invariant models on which we focus in the present 
and the next section K s = 1. For repulsive interactions < K c < 1 while K c > 1 in 
the attractive case. Here we exclusively consider the former. The field operator ^ a (x) 
annihilating an electron with spin direction o =f, \, at position x is decomposed into a 
right- and a left-moving part 

^(x) = e ik ^R a {x) + e- ik ^L a (x). (2) 

The imaginary time fields R a and L a are bosonized according to 

= -y== exp^0 c (r,x)j exp^/ ff ^ B (r,x)J, (3) 



L t(r,x) = -j= exp(^-i0 c (r,x)j exp^-|/ CT s (r, x)J , (4) 



'2-k 

where the Klein factors r] a satisfy anticommutation rules {t]a,Va} = ^o-.o-' and f^ = l = 
—fl- The fields (f> u and (f> u are the chiral components of $^ and Q v , 

®is = <j>v + 4>v, ®v = 4>u - 4>v, v = c,s. (5) 

For the translational invariant TL model the Hamiltonian follows by integrating the 
density Eq. (TT]) over R. The TL model with an open boundary is obtained by integrating 
the density Eq. (|T]) over x > and employing the boundary condition ^f a (x = 0) = 
for the fermionic and $ c>s (x = 0) = for the bosonic fields, respectively. 

We are here interested in the imaginary time-ordered single-particle Green function 

G CTCT /(r,xi,x 2 ) = - (r T ^ C r(r,x 1 )^l / (0,x 2 )) , (6) 
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where (...) denotes the expectation value in the canonical ensemble. From the 
decomposition Eq. ([2]) it follows that 



G n 



D ik F (xi—X2) /~<RR 



_|_ g-i fc F(zi-:r2) QLL _|_ e ik F (x 1 +x 2 ) qRL _j_ e ~ik F (x 1 +x 2 ) Q L R 



where, for example G^£, = — (l~ T i2 CT (r, Xi) ^v(0, x 2 )j. As we are aiming at the local 
spectral function, we eventually set X\ = x 2 = x. For the translational invariant TL 
model the left- and right-moving fermion fields are independent and thus G^, = G^i = 
0. In the presence of an open boundary the above mentioned boundary conditions imply 
R a (r,x) = —L a (r, —x). In this case and after setting x\ = x 2 = x the cross terms are 
characterized by a fast spatial oscillation with frequency 2kp. 

Using standard methods (see e.g. Ref. [2]) one obtains for the Green function of a 
translational invariant system (K s = 1, r > 0, (5 = 1/T, r = x\ — x 2 , R = {x\ + x 2 )/2) 



G™,(t,x 1 ,x 2 ) 



8 n 



271 



7T 



a+b 



7T 



1/2 



Sill 



1/2 



v B /3 



[v.T - ir 



(8) 



X 



sin 

RRi 



[v c T — ir 



sin 



[v c r + ir J 



G%{t, Xi, x 2 ) = G-'{t, x 2 , x h 
and for the case with boundary 



[G^] b (t,x 1 ,x 2 ) 



sinh {j^fjjXi \ sinh ( 



2n 

Vc/3 



X 2 



sm 



[G%] b (t,xi,z 2 )= [G™] b (t,x 2 , Xi ) 



^-Av c t-2iR) sin -^( Vc r + 2iR) 



(9) 



(10) 



The cross terms [G^j] B and [C^] B are equal to — [G^] B and — [G^,] B after 
interchanging r -H- 2R. The appearing exponents are given by 

2 1 / i \ 2 



-4=) , &=I(v^c--t=] , c=!li-^|. (12) 



'#c/ 8 V V^c 

The main steps to obtain the local spectral function from the imaginary time Green 
function are the analytic continuation r — > it + 5 followed by Fourier transformation 
with respect to t. Mathematically the real part 5 corresponds to the ultraviolet cutoff 
introduced to regularize momentum intergrals. From an experimental perspective it can 
be considered as the resolution of the setup (at T = 0). 

For the translational invariant model one obtains ( X \ = x 2 = x) 



p(u),x) 



-(i 



/oo 
. - -oo 



T,X,X) + G^,(r,x,x 



T-Mt+<5 



7T 



2 (1 + 



7T 



2+fe 



7T 



1/2 



sm 



a+fe+l/2 ( 7TT 

/3 



T-Mt+<5 



4 a+fe+1/2 



(13) 
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with 



S y (u) = 2 1+7 r(- 7 ) sin [vr(l + 7)] cosh Q 



1+7 u 



(14) 



« = « + 6-5 = i( Jfc+ i" 2 )- (15) 
The position independent local spectral weight of the translational invariant TL model 
thus shows scaling behavior. T~ a p(cj, x) is a function of uj/T only. The amplitude of the 
scaling function depends on v c / s and K c , which in turn are functions of the interaction 
strength, while its shape is given by K c only. This result was first derived in Ref. [TT] . 
A similar expression was later used to describe transport properties of LLs [27]. 

Taking the T — > limit of Eq. ( TT3]) one obtains the well known power-law 
suppression of the spectral weight [5] EJ [7] 

T = 0: p(u,x) ~ M a , (16) 

for I a; I — >■ 0. For fixed small T > this is cut off by temperature and p saturates for 
\u\ J$ T. For the energy set to the chemical potential, that is u = 0, one finds a power- 
law suppression of p for T — >• (see Eq. (JT3l) ): p ~ T a . From studies of microscopic 
models it is known that (see e.g. Ref. [3]) 

K c = l 




;i7) 



where [/ is a measure of the two-particle interaction and U m is a scale which depends 
on the other model parameters. Using Eq. ( TT5l) one thus obtains 

a ~ f/ 2 , (18) 

for the exponent characterizing the low-energy behavior of the local spectral function 
of the translational invariant TL model. 

Verifying the power-law suppression of the spectral weight as a function of T and 
u with the same exponent a as well as the scaling property of measured STS and/or 
PES data provide strong indications that the system under investigation is indeed a LL 
[22] . It was very recently argued that these characteristics are still not unique to LLs 
as other mechanisms than Id electronic correlations might lead to similar behavior [28J . 
We therefore suggest further consistency checks by in addition measuring spectra close 
to the end points of cut Id chains. 

The local spectral function becomes position x dependent when considering a chain 
with an open boundary. In this case p(u, x) has three contributions 

Pb(u,x) = po(u,x) + e 2ifcF > 2fcp (w,x) + e- 2ifcFX p_2fc F (w,x), (19) 

where the first follows from Fourier transforming G RR + G LL and the last two from 
transforming G and G LR , respectively (x\ = X2 = x). Following the same steps as 
above we obtain 



Po{uj,x) 



-ia+b- 



7T 2 V?+ b vl /2 



- 1/2 F(lu/T,xT/v c ), (20) 
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- y ^ S s in a+b+1 / 2 (i7rs + 8) Vsin(ivrs - 2ti\v + 5) sin(ivrs + 2vrk; + 5)/ 



and 



where 



a+b+l/2 

p 2kF (^x) = p*_ 2fc >,x) = r^-^Giu/^xT/vc), (22) 

2ttV+ 6 ?; s / 



DC 



;i + e -") 1 1 

G(u, v) = as- 



oo 

2c 



(23) 



sin 1/2 (i7rs - 2mvv c /v a + 5) sin a (i7rs - 27rk> + 5) sin fe (i7rs + 2-niv + 5) 
sinh(27ru) 
sin(i7rs + 5) 

In the limit T — > these expressions simplify to the ones given in Refs. [16] and [29] . 

For distances from the boundary beyond the thermal length ~ 1/T, that is 
xT/v c = u > 1 we expect Pb(w, x) to become equal to p(cu, x) of Eq. (fTBl) (exponentially 
fast). That this is indeed the case follows from 
xT 

— > 1 : F(u/T, xT/v c ) « 4S a (u/T), G(w/T, xT/w c ) « 0, (24) 

Eqs. (PI), (120]h fl2"2]) and (jlij). We thus end up with 
xT 

— > 1 : pb(w,x) & p(u,x). (25) 

We next consider the limit xT/v c =ti<l, that is the local spectral function close 
to the open boundary. Then 

/oo ius 
.oo sm a+6+ " !c+i/ ^(i7rs + 5) 
Interestingly the remaining integral has the same form as the one appearing in the 
second line of Eq. (|T3|) but with a + b + l/2 replaced by a + b + 2c + 1/2. For fixed x 
close to the boundary po(x,u) thus displays scaling with the same scaling function as 
the one found in the bulk but a replaced by 

a B = a + b + 2c - 1/2 = i (-L - lj . (27) 

Explicitely one obtains 

— « 1 : po{x,u>) ~ x 2c T a »S aB (u/T). (28) 

With Eq. ( 1TTI) the boundary exponent «b Eq. ( 1271) has, in contrast to the bulk one a, 
a contribution linear in the interaction 

a B ~ U (29) 

and one finds Ob > a. To show that for fixed x close to the boundary T~ aB p-Q indeed 
follows the same scaling function as in the bulk (with a -)• «g) we still have to analyze 
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p2k F for xT/v c = d<1. In the simplest approximation we neglect the v -dependence in 
the integral Eq. (1231) and obtain 



xT 

< 1 : 

Vr 



2Po(^,x). 



(30) 



Using Eq. ( !T9|) this completes our proof that for fixed 

x in the bulk: T' a p(u, x) = T- a p B (co, x) ~ S a (tu/T), (31) 

x close to the boundary: T~ aB p B (u, x) ~ S aB (u/T), (32) 

with Sj given in Eq. (tHJ). 

Showing the consistency of the scaling of spectra measured in the two spatial 
regimes is within reach of the latest STS experiments [22] . Combined with a consistency 
check of the two exponents a and «b, which both depend on K c only and which 
was already achieved in Ref. [22] (see also Sect. @|, this would provide a stringent 
experimental verification of LL physics. 

One can improve the analysis of P2&E close to the boundary by keeping the phase 
factor exp(iftwf) of the integral Eq. (!23l)J|l A numerical evaluation of the integral shows 
that k is a function of K c and v c /v B with n(K c = l,v c /v s = 1) = 2. Taking all terms 
together we find 

1 + «B 



xT 

< 1 : 



T 



p B {u,x) 



Ax 2c cosh (—) 
\2TJ 



■ u \ 



2vrT 



x 



1 



cos 



KUJ 

2k F + — 

Vr 



r 



(33) 



where the overall amplitude A depends on K c , v c and v s but not on the variables u, x, 
and T. 

In the T — > limit Eqs. (fl9"]) to ( 123]) give close to the boundary 



T = 0, 



x\u\ 



< 1 : p B (u,x) ~ x 2c \cu\ aB {l - cos 



2k F + 



KU 



X 



(34) 



At fixed x, pb thus vanishes ~ |w| aB . As in the translational invariant case this power 
law is cut off by a finite temperature and p B saturates for |cu| ~ T. For fixed x close to 
the boundary and uj = Eq. (|32|) gives p B ~ T aB for T — » 0. At T = and deep in 
the bulk, that is for x\u\/v c ^> 1, p±2k F can be written as a sum of terms which vanish 
algebraically in x [TH [29]. They show a power-law dependence on |w| in general each 
with a different exponent. The contribution p to p B becomes position independent and 
goes as po ~ (instead of a B close to the boundary). For sufficiently large x\u\/v c 
(such that algebraically decaying terms can be neglected) one thus finds {uj — > 0) 
x\uj\ 



T = 0, 



> 1 : Pb(w, x) ~ |w| 



For T = and in the noninteracting limit p± 2 fc F does not decay and one obtains 



: 



p B {u,x) 



TTV-p 



COS 



2[kv + — 

V F 



X 



(35) 



(36) 



| We have verified numerically that the additional dependences of the integral Eq. (f23| on u, v, and 
the model parameters are irrelevant for the regimes studied here. 
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for all x and u. 

It is often argued, that the contribution p 2 fc F to pb can be neglected when comparing 
to experiments. The electrons in PES and STS do not come from a specific location 
x but rather from an extended spatial range. If this is large enough compared to the 
characteristic length p2k F averages out due to the fast spatial oscillations with 

frequency 2kp. It is not obvious that the criterion for neglecting p 2fcp is fulfilled in the 
latest STS experiments [22] (see Sect.HJ). 

Equation ( 1331) (or Eq. ( 1341) for T = 0) allows for another consistency check of 
LL behavior. It predicts a spatial power-law dependence of pb close to the boundary 
with exponent 2c = (1/K C — K c )/4 superimposed by oscillations. If it is possible to 
measure the envelope function of the spatial dependence of the spectral weight close to 
a boundary and extract the power-law exponent it would allow to relate the resulting 
K c to the ones obtained from a and/or «b- 

We note in passing that performing a spatial Fourier transform of pb(u,x) reveals 
characteristic informations of the bulk state of a LL including its elementary excitations 
(see Ref. [26] and references therein). 

In the next section we show that the spectral function of microscopic lattice models 
of interacting Id electrons, in our case the extended Hubbard model, indeed shows 
scaling behavior as a function of oj/T. Up to a subtlety in the spatial dependence, 
namely an interaction dependent phase shift in the oscillatory factor, the lattice spectral 
function falls on top of the above computed scaling function of the TL model. This holds 
in the low-energy regime. We discuss the crossover between this universal behavior and 
the nonuniversal regime at higher energies. The crossover scale A depends on the 
parameters of the microscopic model. 

3. Spectral properties of the extended Hubbard model 

Partly aiming at the low-energy LL physics of microscopic lattice models of interacting 
electrons different groups computed the local spectral function [3Ql [31] as well as the 
momentum resolved one [321 [33] for such models using very accurate numerical methods 
(often denoted as 'numerically exact'). Quantum Monte-Carlo (QMC) [301 E3] an d 
(dynamical) density-matrix renormalization group (DMRG) [321 l3Tj were used. The 
results obtained by these methods turned out to be very useful for understanding the 
spectral features of the studied models over the entire band width. Unfortunately, due 
to system size restrictions (DMRG and QMC), artificial broadenings of the spectra 
(DMRG), as well as the problem of analytic continuation of numerical data (QMC), 
it was not possible to reach the low-energy regime; in none of the calculations it was 
possible to convincingly demonstrate power-law behavior of the spectra. Rephrasing 
this in experimental terms one can say that the energy resolution of these methods is 
not high enough|§| For (partly) technical reasons the focus of the numerical approaches 



§ For a very recent attempt to observe power-law behavior of the spectral function in a spinless lattice 
model using DMRG, see Ref. [34]. 
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lies on the Hubbard model with a local two-particle interaction. As the crossover scale 
A between LL behavior and nonuniversal physics in this model is very small (see below) 
reaching the LL regime is particularly challenging. 

We here use a method which allows to obtain approximate results for the spectral 
function of the extended Hubbard model with a local U and nearest-neighbor V 
interaction. Our approximation is based on the functional RG approach to quantum 
many-body physics [25]. Functional RG allows to set up a hierarchy of approximation 
schemes with the two-particle interaction being the small parameter. The one we 
are using here is controlled to leading order and can thus only be used for small to 
intermediate U and V (compared to the band width). The method was mainly applied 
in the context of transport through inhomogeneous LLs and there it was shown to 
reproduce typical impurity strength independent [H] LL exponents to leading order in 
the interaction [35], [36] . Due to appropriate resummations of classes of diagrams the RG 
procedure thus goes way beyond standard perturbation theory. As the exponent of the 
bulk local spectral function is of second order in the interaction (see Eq. (118jl ) the LL 
physics of translationally invariant systems cannot be assessed in our approximation. 
In the following we therefore restrict ourselves to the model with open boundaries 
characterized by the exponent «b which has a linear contribution (see Eq. (|29|) ) and 
study pb close to one of the boundaries. We here refrain from giving any further 
technical details on how the spectral function can be computed within our functional 
RG approach. Those can be found in Ref. [36] . 

The Hamiltonian of the extended Hubbard model with two open boundaries is given 

by 

N-l N N-l 

H = -t Y Y ( c ]+vs> + 4> c j+v ) + u Y n ^ n ^ + v Y n i n J+ i ' ( 3? ) 

3=1 a j=l j=l 

where c\ and c ■ are creation and annihilation operators for fermions with spin a on 
lattice site j, while n Jj(J = cj c-^, and n 3 - = rij^ + is the local density operator 
on site j. For the (nonextended) Hubbard model the nearest neighbor interaction V 
vanishes. The number of lattice sites is denoted by N. The noninteracting tight-binding 
part gives the standard dispersion = — It cos k with the hopping matrix element t > 
(the lattice constant is chosen to be unity). 

Under the assumption that a given microscopic model is a LL (at low energy 
scales) one can use general relations between the exact ground state energy Eq and 
K c [T2l [1} El [3] to extract the dependence of the LL parameter K c on the parameters 
of the model considered. In general however Eq of a many-body system is not known 
analytically. The translational invariant Hubbard model constitutes one of the rare 
exceptions and closed expressions for Eq in form of integral equations can be determined 
using Bethe ansatz [37]. The integral equations can easily be solved numerically which 
gives access to the dependence of K c on the parameters U/t and the band filling n (n can 
vary between and 2) [38] • We emphasize that this only implies that on asymptotically 
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small scales one can expect power-law behavior with a (bulk) or (close to the 
boundary) while no information on the crossover scale A from nonuniversal to universal 
LL behavior can be extracted this way. Furthermore, this expectation holds under the 
assumption that the Hubbard model is a LL, which away from half-filling n = 1— 
for which it is a Mott insulator (37] — is not doubted seriously, but also not proven 
rigorously For the (translational invariant) extended Hubbard model K c (U/t,V/t,n) 
can only be computed numerically along similar lines [39]. The exponents a and Ob 
for the Hubbard [38] cannot become large enough to match the exponents inferred from 
experiments on different systems [18j [191 EH I2H [22] , or, putting it differently, K c cannot 
become small enough. For the extended Hubbard model K^s of roughly the correct 
order can be achieved for U and V of the order of the band width or larger. This part 
of the parameter space lies very close to the Mott transition of the model [39] . One can 
expect that this effects the spectral properties. When aiming at a typical LL spectral 
function with a and «b of experimental size it is thus advisable to study models with 
interaction of longer spatial range. We note that within our approximate approach U 
and V are bound to be sufficiently smaller than the band width. The extended Hubbard 
model is spin-rotational invariant which implies K s = 1. 

For the Hubbard model with an open boundary the scale A was earlier computed 
in the small U limit and it was estimated to be exponentially small [ID] 



In fact, approaching the chemical potential uj — > the local spectral weight first increases 
before the LL power-law suppression sets in for \oj\ J$ A (see the solid line in Fig.[T]which 
does not look 'LL-like' around u = 0; the power- law suppression is beyond the energy 
resolution). This is consistent with the observation of a small crossover scale and a 
peak close to u = in the local spectra of the translational invariant Hubbard model 
obtained numerically by QMC and DMRG [501 151] . 

To compute the finite temperature p-B(oo,j) (here the continuous position x is 
replaced by the discrete lattice site index j) of the extended Hubbard model we consider 
a chain of N lattice sites described by Eq. (15T1) . For this the spectrum is discrete and 
the spectral function consists of <5-peaks. Due to even-odd effects the spectral weight 
might vary quickly from one eigenvalue to the next one. A smooth function of u is 
obtained by averaging the weight over neighboring eigenvalues}]]] To obtain the local 
spectral function as defined in the continuum one furthermore has to devide the weights 
by the level spacing between eigenvalues. The energy scale 5n = nvp/N associated to 
the chain length becomes irrelevant as we always consider sufficiently large systems with 
T ?t 5n for fixed T. Our results are thus not influenced by finite size effects (for an 
exception, see the discussion of Fig. Hj). Typical experimental temperatures are in the 
few to few ten Kelvin range which corresponds to T « 10~ 4 t to 10~H for our model. 

|| A similar energy averaging is inherent to any STS or PES experiment due to the finite energy 
resolution of these techniques. 




(38) 




qU , 1 , 1 , 1 1 

-1 1 2 

co/t 

Figure 1. (Color online) Local spectral function of the extended Hubbard model for 
n = 3/4, j = 1, U/t = 0.5, N = 2 14 , T/t = 10~ 3 and different V/t. For filling n = 3/4 
the optimal nearest-neighbor interaction is given by V /t = U/(tV2) « 0.35. Only for 
sizable nearest-neighbour interaction V we observe the LL suppression of the weight 
at uj = 0. The suppression close to u/t = 1.5 is a lattice effect [4Tj . 



In Fig. [T] pb(^,j) is shown for filling n = 3/4, lattice site j = 1 next to the 
boundary, U/t = 0.5, iV = 2 14 , T/t = 10 -3 and different V/t. For the Hubbard model 
with V/t = no suppression of the spectral weight is observable as A is much smaller 
than temperature. Obviously, A increases with increasing V/t and the LL suppression 
at uj = becomes apparent. This can be understood as follows. The crossover scale 
A is strongly affected by the size of the open boundary analog of a g± : ± two-particle 
scattering process [T71 [361 HO] which cannot be written quadratically in the bosonic 
densities. Its initial value (with respect to an RG flow) in the extended Hubbard model 
is given by g± : ± = U + 2V cos(2/;; F ) with kp = mi/2. It is large for V = 0. At 
V = — U/[2 cos(2/cf)] it vanishes (see the dashed-dotted line in Fig. []]). Under an RG 
procedure this 'non-LL term' flows to zero and is thus RG irrelevant. However, the flow 
is only logarithmically. This implies that for sizable initial g± ± LL physics sets in on 
exponentially small scales consistent with Eq. ( 1381) for the Hubbard model [10]. For 
small initial <7i,_i_, Pb appears LL-like with the characteristic power- law behavior of the 
spectral weight close to u = 0. In Fig. [T] the spectral weight at uj = remains finite due 
to the finite temperature. We conclude that to observe LL physics on moderate scales in 
the present model the interaction should not be too local. In particular, to demonstrate 
power-law behavior and obtain an estimate of the exponent by fitting pe for fixed j 
(close to the boundary) and T C f as a function of u (in the range T \cu\ Ct) one 
should consider fine-tuned parameters with V = V a . The spectral functions of Fig. CD 
show another 'high-energy' nonanalyticity. A similar feature was observed for a model 
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Figure 2. (Color online) Low-energy regime of the local spectral function of the 
extended Hubbard model for n = 3/4, U/t = 0.5, V = V a = U/y/2 (for n = 3/4), 
_ 2 14 . T/t = 10~ 3 and different j close to the open boundary at j = 1. 



of spinless fermions in Ref. [JT] and was explained there as a lattice effect. 

In Fig. |2] the low-energy regime of pe is shown for the same parameters as in Fig. [1] 
but 'optimal' V = V Q = U/\/2 (for n = 3/4), which allows for the largest low-energy 
regime, and varying position j close to the boundary site j '• = 1. For fixed u we observe 
strong variations of the weight with j and pronounced u -H- — w asymmetries which is 
consistent with the result from the TL model Eq. (I33p . Below we return to the spatial 
dependence of p-Q. 

To confirm scaling in oj/T at fixed j as predicted in Eq. fl3"3"|) we computed p& for 
the parameters of Fig. |2] but with j = 1 and for different T. By fitting p^(oj,j = 1) 
as a function of u for the smallest T in the range T \uj\ £ we can extract a 
functional RG estimate of the boundary exponent a^ G . For the given parameters we 
obtain a^ G = 0.089 in good agreement with the DMRG result ag MRG = 0.095 obtained 
from Eq. ( |27j) and K® MRG = 0.840 derived as explained above [39]. The scaling obtained 
with this is shown in Fig. [3] In the inset the unsealed data are displayed. The thick 
solid line is the prediction Eq. ( |33|) of the TL model with an open boundary, where we 
replaced Ob —> a^ G . The data nicely collapse on the TL model curve in the low-energy 
regime. 

We next take a closer look at the j dependence of pe at fixed u. As emphasized 
in the last section measuring the local spectral weight as a function of j offers another 
possibility for a consistency check that the system under consideration is a LL: Eq. (1331) 
predicts power-law behavior of the envelope with exponent 2c = (1/A' C — K c )/4. Figure 
S shows p B for n = 3/4, U/t = 0.5, N = 2 14 , T/t = 10" 3 and u w as a function 
of j (filled circles). Here oj ~ refers to taking the eigenvalue of the finite system 
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Figure 3. (Color online) Scaling plot of the local spectral function of the extended 
Hubbard model for n = 3/4, U/t = 0.5, V = V a = U/V2 (for n = 3/4), TV = 2 14 , j = 1 
and different T. The inset shows the unsealed spectral function as a function of ui for 
different T. 



closest to u — 0, which might be of order 1/tV away from zero. We again tune 
V to the optimal value V Q = U/y/2 (for n = 3/4) providing the largest A. The 
spatial oscillations with frequency 2/c F = 37r/4 are apparent. We fitted the envelope 
to a power law ~ j^c mG an( j obtained 2c fRG = 0.087 in excellent agreement with 
2c dmrg = (i/j^DMRG _ k? mrg )/4 = 0.088. The power-law fit is shown as the thick 
solid line in Fig. HI As mentioned above we control the different exponents only to 
leading order in the interaction. To this order the analytic expressions for ob and 2c 
agree, as is apparent from Eqs. (I12p . (T27j) and (fT7|) . The numerical values for 2c and 
a^ G still differ by roughly 2% as the RG produces higher than linear order terms in 
the different exponents as well. We emphasize that for V — 0, that is for the Hubbard 
model, in a similar plot no spatial suppression of the envelope of the spectral weight at 
small j is visible. In fact, the envelope of the spectral weight at u ~ even increases 
for j approaching the boundary site j — 1. 

A comparison of the data for U/t = 0.5 (filled circles) and for U/t = (crosses) 
additionally presented in Fig. |4] shows that a phase shifty of the spatial 2kp oscillations 
appears which is not captured by the result for the TL model Eq. (J33l) . The latter was 
derived using standard bosonization with the boundary conditions on the continuum 
fields given after Eq. (jSJ). The deviation of the U/t = curve from Eq. (1361) for j <J 20 
(splitting of degenerate values of the spectral weight) is a finite size effect; u is not 
exactly zero but of order 1/tV (eigenvalue closest to zero). The phase shift can be most 
easily identified from the observation that pe vanishes on every eighth lattice site for 
U/t = (crosses), as it is supposed to according to Eq. (1361) with oj = and 2/c F = 37r/4, 
but not so for U/t — 0.5 (filled circles), where the same should hold according to Eq. (|33|) . 
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Figure 4. (Color online) Lattice site dependence of the spectral function of the 
extended Hubbard model close to the boundary site j = 1 for n = 3/4, two different 
U/t, V = V = U/y/2 (for n = 3/4), N = 2 14 , T/t = 1(T 3 and u w 0. The solid line is 
a power-law fit to the envelope. 



The phase £ turns out to be linearly dependent on U (for small U/t and V = V Q {U)). 
As our functional RG approximation scheme is controlled to this order, the appearance 
of £ is a reliable finding. Considering T = at different system sizes N we furthermore 
verified that £ does not vanish for decreasing 1/N. The phase shift is thus not a finite 
size effect. A phase shift as observed in the extended Hubbard model (with V = V (U)) 
can be accounted for in bosonization by adding a local single-particle forward scattering 
term WS(x) d x Q c (x) to the Hamiltonian density Eq. (pQ). The phase shift is then given 
by 2nK c W/v c . To match the result of the extended Hubbard model (with V = V Q (U)) 
W has to be chosen [/-dependent, in particular W ~ U for small U/t. The STS and 
PES experiments always imply a spatial averaging. As the phase shift becomes illusive 
even after averaging over only a few lattice sites and as the main focus of the present 
paper is on relating theoretical spectral functions to experimental ones we here do not 
further investigate this issue. 

Finally we compare the spatial dependence of pb for two different u in Fig. [5j 
Apparently the frequency of the spatial oscillations depends on u which is consistent 
with Eq. (133|) . As this is more transparent for a 'more commensurable' filling, the 
parameters of this figure are n = 2/3, U/t = 0.5, V — V — U (for n = 2/3), N = 2 14 
and T/t = 1(T 3 . 
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Figure 5. (Color online) Lattice site dependence of the spectral function of the 
extended Hubbard model close to the boundary site j = 1 for n = 2/3, two different 
w, U/t = 0.5, V = V = U (for n = 2/3), N = 2 14 , T/t = 10~ 3 . 

4. Summary 

In the first part of this paper we have derived analytic expressions for the power- 
law behavior of the local spectral weight p of the translational invariant TL model 
as well as the one with an open boundary p^ as a function of u, T and x. The results 
provide a variety of possibilities for consistency checks of experimental STS and PES 
data on Id electron systems. The first is to show scaling of data for different u and 
T taken in the bulk or at fixed position close to the boundary onto the bosonization 
predictions Eqs. (fl3|) and (TL31) (bulk) or Eq. (|33|) (boundary). Scaling of bulk spectra 
was e.g. demonstrated in Ref. [22]. If the same could be achieved for boundary spectra 
the same scaling function with a replaced by «b should appear if the system is a 
LL. Experimentally showing this together with the required consistency of a and Ob 
(for a spin-rotational invariant model with K 8 — 1 both given by a single number 
K c ; see Eqs. (I15p and (1271) ) would constitute a second highly nontrivial check that the 
studied system indeed is a LL. A consistency of bulk and boundary exponents within the 
experimental error bars (but not the entire scaling function) was achieved in Refs. [27J 
and [22] (also see below). A third consistency check is provided by the predicted spatial 
power-law behavior with the exponent 2c which again can be expressed solely in terms 
of K c (see Eq. (112])). 

It is often argued that the spatially oscillating contributions pik ¥ and p-2k F to pe 
with frequency 2kp can be neglected when comparing to experimental spectra due to 
spatial averaging effects. Taking the numbers from Ref. [22] this is not apparent. In 
this experiment hp ~ 5 x 10 8 /m while the range of spatial averaging is estimated as 
Ax ~ 5 x 10 _9 m, leading to 2kpAx = 5. This is not a very large number and one 
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would thus conclude that p 2 fc F and p-2k F cannot be dropped. A quantitative picture of 
the averaging effects can easily be obtained by integrating Eq. ffT9~]) over an appropriate 
spatial range. We have explicitely verified that averaging over Ax ~ 5 x 10~ 9 m does not 
significantly smear out the boundary (xu/v c <C 1) and bulk (xu/v c ^> 1) exponents in 
oj at low temperatures (here T = 0). In particular, this shows that the spatial resolution 
of the experiment is high enough to detect a B (as it is implicit to the analysis presented 
in Ref. [22]). 

In connection with the comparison of the experimental results on gold chains on 
a germanium surface of Ref. [22] to the LL predictions one might be worried about 
two effects which are not included in the TL model of Sect. [2j One is Rashba spin- 
orbit interaction (SOI), which in a surface setup can become sizable. Along the lines of 
Refs. [42] and [43] one can bosonize the Id electron gas with Rashba SOI. An important 
effect is the appearance of two different Fermi velocities up(l±£) due to subband mixing 
and the SOI splitting j32J EE]. Here ( is a measure of the strength of the SOI. The local 
spectral function (with and without an open boundary) shows the same characteristics 
as a function of to as discussed in Sect. [21 but with modified exponents 



aS0 ' — W^t) ( + 0{<% (39) 

«I 0, = «b + ^C + o(C 2 ). (40) 

Taking realistic numbers for the velocities it turns out that ( < 1 and the effects of 
Rashba SOI on the exponents are negligible. We note in passing that also the momentum 
resolved spectral function of a translational invariant LL is barely modified by SOI of 
realistic size |44j . 

The other issue is the observation of four electron branches (instead of two in the 
TL model) crossing the Fermi surface in PES measurements on the gold chains [35]. The 
four branch situation can also be accounted for in bosonization [4"6| 147] . One then has 
to introduce even and odd pairs of charge and spin density bosons as well as the related 
K's and velocities. Under the plausible assumption that only the even charge channel 
LL parameter is different from the noninteracting value 1, that is K ce — K c < 1, one 
finds the same power-law behavior for the local spectral functions as a function of u as 
the ones given in Sect. [2] but with 

a =K^ + i^ 2 )' (41) 

<*B = 7 ( TT - 1 ] ■ (42) 



,2 



4 V^c 

Interestingly taking these expressions would significantly improve the consistency of the 
experimentally determined a and «b [22] ■ From the measured value a = 0.53 of the 
bulk spectra one finds Ob = 1.27 which nicely agrees to the measured value a-Q = 1.20. 
In particular, the agreement is improved compared to taking the two-branch expressions 
Eqs. (EES]) and (EZJ) which gives a B = 1.43 [22]. 
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In the second part of our paper we have shown that the behavior of pb as a function 
of u, T and j close to an open boundary predicted by the bosonization solution of the 
TL model can indeed be found in an example of a microscopic lattice model, namely 
the extended Hubbard model. Interestingly, a linear-in-[/ phase shift £ not captured by 
standard open-boundary bosonization appears in the spatial oscillations of the spectral 
weight of the extended Hubbard model. The energy resolution required to access the 
low-energy LL regime strongly depends on the model parameters. Even for the fairly 
high energy resolution we can achieve within our approximate method, convincingly 
demonstrating LL power-law behavior and reliably extracting exponents requires fine- 
tuning of the parameters. Roughly speaking the crossover scale becomes small if the 
two-particle interaction becomes too local. This suggests that to access the low-energy 
LL regime at a given experimental energy resolution one should consider systems with 
poor screening properties. 
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